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PROCEDE DE RESTAURATION SPATIALE DES COULEURS D'UNE 

IMAGE • 

La presente invention appartient au domaine technologique de 
l'imagerie. La presente invention a pour objet un procede autoraatique de 
restauration de la couleur d'une image, apres que ladite couleur ait ete alteree par 
rapport a la couleur d'origine de l'image. 

La renovation ou restauration des images en couleur enregistrees sur 
des supports argentiques, est generalement operee en faisant un traitement d'une 
version numerisee de ces images argentiques alterees. La restauration des images 
alterees permet qu'a l'examen visuel, l'image recouvre ses couleurs d'origine. Les 
supports traditionnels argentiques des images sent constitues generalement par du 
film photographique, ou du papier photographique. On designera dans la 
description qui suit, par image argentique, une image en couleur obtenue a partir 
d'emulsions photographiques aux halogenures d'argent. L'alteration ou degradation 
des images en couleur enregistrees sur ces supports argentiques est due par 
exemple au vieillissement, ou a 1'exposition de ces supports argentiques a des 
sources lumineuses, a l'occasion de leur manipulation. Les procedes de restauration 
des couleurs d'une image, connus de Homme du metier, consistent par exemple a 
numeriser une image prealablement emegistree sur un support argentique (film ou 
papier photographique), puis a traiter l'image numerisee afin de restaurer ses 
couleurs. Ces procedes de restauration consistent a transformer les couleurs 
degradees sur l'image numerisee, en mettant en oeuvre des algorilhmes de logiciels 
de traitement d'image. 

L'image numerique, transformee par le proced6 de traitement, 
presente ainsi des couleurs restaurees qui correspondent approximativement a 
celles de l'image d'origine. 

Le brevet US 5 796 874 decrit la restauration d'images imprimees 
sur un support materiel, par exemple de type papier, decolor6es ou fanees au cours 
du temps. Ce brevet decrit une methode et des moyens qui permettent a un 
operateur n'ayant pas un niveau de competence eleve, de restaurer et facilement 
une image dont les couleurs ont ete alter6es. L'operateur n'a notamment pas a 
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choisir et a ajouter des filtres dans le scanner utilise pour num6riser l'image a 
restaurer. Mais, 1'operateur doit appliquer des modeles de degradation de l"image 
au cours du temps a l'image a restaurer. C'est-a-dire que la restauration est operee 
en fonction du choix d'une variable de temps: par exemple, on suppose que l'image 

5 a une dizaine d'annees d'age et on applique le modele de restauration correspondant 
a ce temps (une dizaine d'annees) afin de restaurer les couleurs de l'image alteree. 
La methode de restauration permet d'obtenir et de selectionner automatiquement un 
modele de restauration des couleurs de l'image, en fonction d'une variable 
temporelle affectee a l'image. L'image ainsi restauree peut 6tre affichee sur un 

10 ecran de type moniteur. De plus, 1'operateur peut aussi choisir une autre modele, lie 
a un autre variable de temps, si le premier resultat de restauration ne le satisfait 
pas. La regie de correction utilisee est appliquee uniformement a l'ensemble de 
l'image. 

Le document de lUniversite de La Rochelle, France, de M. 

15 Chambah et B. Besserer, intitute Digital Color Restoration of Faded Motion 
Pictures, et presente a la premiere Conference Internationale CGIP (Color in 
Graphics and Image Processing) de Saint-Etienne, France, du 01 au 04 octobre 
2000, decrit une methode de restauration des couleurs damages fanees sur de vieux 
films. Dans un premier temps, le film est immense avec un scanner; puis l'effet du 

20 recouvrement des densites spectrales des differents colorants est soustrait, en 

utilisant une matrice d'ajustement Dans un second temps, les canaux de couleurs 
de l'image sont equilibres, en utilisant une autre matrice de correction. Finalement, 
le contraste est augmente, pour ameliorer la qualite visuelle. Dans cette methode, 
la regie de correction utilisee est appliquee uniformement a l'ensemble de l'image. 

25 Dans un article (IEEE Transactions On Image Processing, volume 9, 

numero 5, de mai 2000), intitule "Adaptive Image Contrast Enhancement Using 
Generalizations of Histogram Equalization", J. Alex Stark decrit une methode 
d'augmentation du contraste dans une image. Cette methode utilise une fonction 
cumulative des niveaux de gris dans une zone ou une fen&re (par exemple de 

30 forme carr6e) autour d'un pixel, sans contrainte d'interdependance entre ces zones. 
Les corrections sont independantes, pour chaque zone de representation de l'image. 
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Les moyens de restauration de la technique anterieure ne prennent 
pas en compte les variations spatiales de l'alteration des couleur dans une image. 
Les moyens de la technique anterieure utilisent des modeles ou fonctions 
mathematiques permettant de prendre en compte certains parametres importants 
qui influenced la degradation des couleurs d'une image enregistree sur un film ou 
un papier photographique, pour ensuite restaurer les couleurs de l'ensemble de 
1'image d'une maniere uniforme. 

Le proceae de la presente invention, pour la restauration spatiale des 
couleurs d'une image alteree, permet, a partir de la detection des alterations ou 
degradations de couleur qui varieht spatialement d'un bord a l'autre de 1'image, une 
restauration des couleurs conformement aux couleurs d'origine de ladite image. 
L'avantage de ce procede de invention est de restaurer automatiquement, 
rapidement, et de maniere fiable les couleurs d'origine de 1'image, alors que lesdites 
couleurs ont et6 alterees non uniformement. 

Plus particulierement, l'invention concerne un precede" automatique 
de restauration de la couleur d'une image argentique dont la couleur est alteree 
d'une facon non uniforme par rapport a la couleur d'origine de 1'image argentique, 
spatialement dans le plan de 1'image argentique, et suivant une direction parallele a 
un axe principal de 1'image; ce proced6 comprend les etapes suivantes: 

a) numeriser 1'image argentique alteree en utilisant une echelle de 
mesure des densites optiques; 

b) filtrer l'ensemble des pixels de 1'image numerisee; 

c) diviser 1'image numerique en bandes de pixels disposees 
perpendiculairement a la direction de l'axe principal de 1'image 
suivant lequel se produit l'alteration de couleur, la reunion de 
l'ensemble des pixels de chaque bande representant l'ensemble 
des pixels de 1'image numerique; 

d) calculer. pour chaque bande de pixels, et pour chacun des 
canaux de couleur formant 1'image, une distribution des densites 
optiques des pixels formant ladite bande; 
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e) calculer, pour chacune desdites distributions de l'etape d), une 
valeur de reference brute superieure de densite optique, et une 
valeur de reference brute inferieure de densite optique; 

f) calculer, a partir de l'ensemble des valeurs de reference brutes 
superieures et inferieures de l'etape e), des valeurs de reference 
filtrees superieures de densit6 optique et des valeurs de 
reference filtrees inferieures de densite optique correspondantes; 

g) determiner, a partir des valeurs de rdference filtrees superieures 
et inferieures de density optique de l*6tape f), des courbes de 
reference superieures de densite optique et des courbes de 
reference inferieures de densite optique, lesdites courbes de 
reference etant representatives du profil de degradation de 
couleur de l'image; 

h) determiner, pour chacun des canaux de couleur formant l'image, 
et a partir de la valeur maximale de densite optique de la courbe 
de reference superieure de densite optique, la bande de pixels de 
moindre degradation; 

i) appliquer a chaque bande de pixels de l'image numerique autre 
que la bande de pixels de moindre degradation, et pour chacun 
des canaux de couleur formant l'image, une transformation 
lineaire permettant de transformer les valeurs filtrees de densite 
optique placees sur la courbe de reference superieure et sur la 
courbe de reference inferieure de densit6 optique de ladite 
chaque bande de pixels, respectivement en valeurs superieure et 
inferieure de densite optique de la bande de pixels de moindre 
degradation. 

Le precede de invention permet d'obtenir une image restauree. Une 
caracteristique de Invention qui se traduit par un avantage majeur par rapport a la 
technique anterieure, est de prendre en compte les variations spatiales non 
uniformes de couleur ayant altere par exemple une image argentique. Ce proced6 
permet de restaurer automatiquement les couleurs fan6es de l'image, sans faire 
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intervener les competences photographiques d'un operateur. Un autre avantage du 
precede de Tinvention pour executer des restaurations de couleur d'images alterees 
ou fanees est le temps gagne, et par consequent une productivity accrue, par 
rapport a un traitement manuel de la restauration spatiale des images par un 
operateur. 

D'autres caracteristiques et avantages apparaitront a la lecture de la 
description, faite en reference aux dessins des differentes figures. 

La figure 1 represente schematiquement un film photographique en 
couleur et un moyen destinS a son enroulement 

La figure 2 represente schematiquement un exemple 
d'environnement materiel (hardware) dans lequel est par exemple mis en ceuvre le 
proced6 de Tinvention. 

La figure 3 represente schematiquement une image num6ris6e de 
Timage a restaurer. 

La figure 4 represente un exemple de distribution, dans un canal de 
couleur, des densites optiques d'une bande de pixels d ! une image numerique. 

La figure 5 represente schematiquement, pour un canal de couleur, 
une Evolution de la distribution des densites optiques dans Vimage, selon un axe 
principal suivant lequel se produit Talteration de couleur. 

La figure 6 represente schematiquement un exemple de diagramme, 
dans chaque canal de couleur, de representations spatiales de la distribution des 
densites optiques dans Timage filtree de son contenu. 

La figure 7 represente la fonction de restauration des densites 

optiques. 

La figure S represente un exemple de distribution des densites 
optiques cumulees de la figure 4. 

La figure 9 represente un principe de filtrage de Vimage par suivi 

des maximum. 

La figure 10 represente schematiquement, pour un canal de couleur, 
la distribution des densites optiques apres filtrage, de la figure 5. 



La description suivante est une description detaillee de l'invention 
en reference aux dessins dans lesquels les memes references numeriques identifient 
les memes elements dans chacune des differentes figures. 

Le precede de l'invention est utilise pour restaurer automatiquement 
les couleurs dhine image dont la couleur est alteree, c'est-a-dire fan6e par rapport a 
la couleur d'origine de l'image. Les images a restaurer peuvent etre fixes ou 
animees. Les images fixes sont par exemple des photographies en couleur 
(portraits, paysages, etc.). Les images animees sont par exemple des sequences de 
film cmematographiques. Ces images sont enregistrees sur des supports 
photographiques argentiques. Le probleme d'alteration non homogene des couleurs 
des images enregistrees sur un support photographique argentique est rencontre 
notamment lorsque, suivant la figure 1, ces supports d'images argentiques 2 sqnt 
enroules sur des bobines 1 . Le support d'images argentiques est par exemple un 
film cinematographique enregistre en couleur. Au cours du temps ou suite a des 
variations de temperature et/ou dlrygrometrie non maitrisees, des degradations de 
colorants se produisent. Pratiquement, les degradations des colorants de l'image se 
produisent plus fortement sur les bords 4F, 5F qu'au centre du film 2. La ligne 
materialisant le milieu du film 2 et des images argentiques 3 est representee par un 

premier axe principal 6. 

La figure 2 repr6sente un environnement materiel preferentiel dans 
lequel le precede de l'invention est mis en anrvre de maniere interactive. Un 
operateur qui desire par exemple restaurer une ou plusieurs images d'origine 3 dont 
les couleurs ont ete alterees au cours du temps par rapport a la couleur d'origine, 
utilise un terminal 7. Le terminal 7 est par exemple un ordinateur personnel ou PC 
(Personal Computer) comprenant un 6cran d'affichage 1 1 . Le terminal 7 est 
interface a une souris 9 comprenant au moins un bouton de commande 10. Le 
systeme d'exploitation du terminal 7 est par exemple Unix, Windows, ou Linux. 
Tous ces systemes d'exploitation permettent de reconnaitre des peripheriques 
comme un scanner ou numeriseur 13, une imprimante 14, ou tout autre dispositif 
15 rehe au terminal 7, et permettant la capture d'images ou de sequences d'images. 
Le proc&le de l'invention permet de numeriser automatiquement l'image 




argentique 3 dont les couleurs sont alterees, en utilisant par exemple le scanner 13. 
L'image 3 numerisee est codee en valeurs numeriques. Pour numeriser l'image, on 
utilise une echelle de mesure des densites optiques. Autrement dit, les valeurs 
numeriques associees a chaque pixel, et pour chacun des canaux de couleur 
formant l'image, sont proportionnelles (fonction lineaire) aux densit6s optiques. 
L'image 3 est numerisee en une image 12 qui peut ainsi 6tre visualis6e sur l'ecran 
1 1 . L'image numerique 12 peut etre memorisee dans une memoire interne au 
terminal 7. 

Dans un autre mode de realisation, la mise en ceuvre du precede de 
l'invention peut aussi ne pas etre interactive, et n r utilise pas forc&nent des moyens 
comprenant un ecran d'affichage 1 L Dans ce mode de realisation, le precede de 
Invention est mis en oeuvre en utilisant un processeur dedie, c'est-a-dire un 
processeur configure avec des portes et des fonctions logiques de base permettant 
de mettre en oeuvre ledit precede. On utilise par exemple un processeur FPGA 
(Field Programmable Gate Array). Ce processeur est par exemple integre dans le 
scanner 13, et execute automatiquement la restauration juste apres l'operation de 
numerisation de l'image argentique 3. 

L'alteration, c'est-&-dire la decoloration ou le fanage des couleurs 
d'une image argentique est causee en particulier par l'exposition du support de 
l'image argentique a la lumiere et/ou a MrumiditS. Les films sont par exemple 
stockes enroules sur des bobines 1; dans ces conditions, l'alteration des couleurs de 
1'image argentique 3 enregistree sur le film 2 se produit d'abord au voisinage des 
bords 4F, 5F dudit film 2. Au cours du temps, Talt6ration de couleur est done plus 
importante au voisinage des bords 4F, 5F du film 2 enroule sur la bobine 1, qu'au 
voisinage de l'axe principal 6 situS au centre du film 2. En consequence, une 
premiere hypothese fondee sur des resultats experimentaux, et preaiable a la mise 
en oeuvre du precede de l'invention, consiste k admettre d'une part que l'alt&ration 
de couleur varie spatialement dans le plan de l'image 3, suivant une direction 
perpendiculaire au premier axe principal 6, ou parallele k un second axe principal 
19 de l'image 3, 12. Done, Alteration de couleur varie en fonction de la distance 
par rapport aux bords 4F, 5F du film 2. D'autre part, cette premiere hypothese est 
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liee a une seconde hypothese qui consiste a admettre egalement que la degradation 
de couleur est uniforme (c f est-a-dire ne varie pas) suivant le premier axe principal 
6, parallele aux bords 4F, 5F. Cette hypothese d^niformite de la degradation de 
couleur doit etre verifiee au moins sur l'etendue de l'image selon cet axe principal 
5 6. 

Selon la figure 3, et par analogic avec la figure 1, les bords 16, 17 
de 1'image num6rique 12 (image 3 numerisee) correspondent respectivement aux 
bords 41, 51 de rimage argentique 3. Les bords 26, 27 de rimage numerique 12 
correspondent aux deux autres bords de rimage argentique 3, lesdits autres bords 

10 etant paralleles au second axe principal 1 9. Le premier axe principal 6 est 

perpendiculaire au second axe principal 19. Les deux axes principaux 6, 19 de 
Timage 3, 12 se coupent en un point 20 qui est place au centre de rimage. 

Pour determiner les variations des densites de couleur spatialement 
dans le plan de 1'image numerique 12 par rapport a la couleur d'origine de rimage 

15 argentique 3 correspondante, l'algorithme du precede de rinvention permet de 
calculer automatiquement des distributions de pixels, en fonction des valeurs 
numeriques codees de leurs densites optiques respectives, dans une bande mobile 
18 de rimage 12. Le procede de rinvention permet de calculer ces distributions 
dans des tranches ou bandes 18 representatives d'une partie de rimage 12. La 

20 largeur 21 de la bande 18 inclus par exemple cinq colonnes de pixels paralleles a 
l r axe 6. La longueur de la bande 18 correspond a l'extension totale de rimage 12 
suivant la direction de l'axe 6, perpendiculaire a la direction de l'axe 1 9; la 
direction de l'axe 19 caracterise la direction de la variation de l ! alteration des 
couleurs. ^extension totale de rimage 12 est la dimension de rimage: par exemple 

25 sa hauteur ou sa largeur. Dans un mode de realisation particulier de rinvention, la 
bande 181a moins large est reduite a une seule rangee de pixels paraltele a l ! axe 6. 

La figure 4 repr6sente un exemple de distribution des densites 
optiques d'une bande de pixels, dans un canal de couleur donne, par exemple le 
rouge. Une courbe 22 represente la distribution des pixels (nombre de pixels N 

30 suivant l r axe des ordonnees) en fonction des valeurs codees de densit6 optique 

V c (DO) desdits pixels (suivant l ! axe des abscisses). La courbe 22 represente, pour 
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une position spatiale donn6e de la bande 1 8 dans le plan de l'image 12, la 
distribution des valeurs codees de densite optique Vc(DO) des pixels formant 
ladite bande 18, entre une valeur niinimale de densite optique representee par le 
point 23 et une valeur maximale de densite optique representee par le point 24. 
Generalement, Pimage est affectee par du bruit (signal parasite). Les causes de 
bruit peuvent etre multiples. Dans le contexte de Pinvention, le bruit genant est du 
principalement a des poussieres. Les poussieres ont g£ner6 par exemple des points 
noirs dans Pimage, ou bien des points blancs, et perturbent la distribution des 
valeurs codees de densite optique Vc(DO). Selon la figure 4 9 les valeurs codees de 
densite optique des points noirs sont par exemple representees par un lobe 28, et 
les valeurs codees de densite optique des points blancs sont par exemple 
representees par un lobe 29. Le procede de ^invention permet d'eliminer 
automatiquement ces lobes aberrants 28, 29. 

La figure 8 represente un exemple de courbe 70, de la distribution 
des densites optiques cumulees de la figure 4. Selon un mode de realisation 
particulier de Pinvention, la distribution cumulee de la figure 8 est par exemple 
divisee en classes d'effectifs egaux de pixels ayant des valeurs voisines de densites 
optiques. Preferentiellement, la distribution est divisee en centiles (100 classes). 
Mais la distribution peut etre divisee en un autre nombre arbitraire de classes. Pour 
eliminer les lobes aberrants eventuels 28, 29 dus au bruit, le procede de Pinvention 
ne traite que la partie de la courbe 70 comprise entre les points d'ordonnees 
respectives X x et (100 - X'0. Dans un mode de realisation prefere: Xi = X\ = 5%. 
Mais les valeurs de Xi et X ! i sont parametrables, c'est-a-dire peuvent prendre 
d'autres valeurs (pourcentage), et ne pas §tre egales entre elles. 

Une troisieme hypothese de mise en ceuvre du procede de 
Pinvention consiste a considerer que le changement de densite optique 
correspondant a Talteration de couleur s f opere suivant un modele lineaire, 
caracterise par exemple par une premiere Equation Eq(l) du type: D F = a Dr + p. 

Le mode de realisation prefere de realisation de Pinvention utilise le 
modele lineaire en Pappliquant aux densites optiques de pixels de bandes 1 8 de 
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l'image numerique 12. Dans un autre mode de realisation, on peut egalement 
utiliser ce modele lineaire en l'appliquant aux intensites des couleurs. 

Dans I'equation D F = a Dr + p, D F represente la valeur codee de la 
densite optique d'un pixel altere, c'est-a-dire un pixel dont la valeur codee de 
5 densite optique a par exemple diminue par rapport a la valeur codee de densite 
optique de reference Dr dudit pixel. La valeur codee de densite optique de 
reference Dr correspond a la densite optique de la couleur d'origine du point 
correspondant spatialement dans le plan de l'image argentique 3, au pixel de 
l'image numerique 12. Les valeurs de a et P varient selon l'axe principal 19 de 

10 l'image numerique 12. 

Une quatrieme hypothese prealable a la mise en oeuvre du procede 
de l'invention, admet que le profil de degradation r6el de couleur de l'image varie 
de maniere continue (c'est-a-dire sans rupture), lisse (c'est-a-dire avec une variation 
lente) et monotone (c'est-a-dire avec une degradation toujours decroissante du bord 

15 vers le point de moindre degradation). Le point de moindre degradation de l'image 
12 est represente par exemple par le point 20 place au centre de l'image 12. Mais 
ce point de moindre degradation peut gtre place spatialement n'importe ou dans le 
plan de l'image 12. Le profil de degradation detecte a pour objectif de se 
rapprocher du profil de degradation reel. Done, le resultat de restauration de 

20 couleur est d'autant meilleur que le profil de degradation detecte par le precede de 
l'invention est continu, lisse, et monotone. C'est-a-dire que le profil s'affranchit des 
variations cbromatiques inherentes au contenu mSme de l'image. Des moyens de 
filtrage connus de lliomme du metier, permettent de rendre le profil de degradation 
detecte, continu, lisse, et monotone. 

25 Dans un mode de realisation prefere de l'invention, les etapes de 

detection de Talteration des couleurs de l'image sont generalement precedees d'un 
filtrage initial automatique de 1'ensemble des pixels de l'image numerique 12. Le 
filtrage initial permet de minimiser Tinfluence du contenu de l'image et des 
poussieres (artefacts dans l'image). Mais, un autre mode de realisation particulier 

30 permet de ne pas faire ce filtrage initial. Le filtrage de l'image numerique 12 est 

opere en utilisant des operateurs de morphologic mathematique connus de lTiomme 
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du metier. Pour chaque pixel de Vintage 12, est operee une ouverture (erosion 
suivie d f une dilatation) puis une fermeture (dilatation suivie dhine erosion), en 
utilisant un element structurant ayant par exemple une forme carree et contenant 9 
ou 25 pixels. Ces ouverture et fermeture peuvent etre operees consecutivement n 
5 fois, n etant un nombre entier. 

Le procede de rinvention permet de diviser automatiquement 
Timage numerique 12 en bandes de pixels 1 8 disposees perpendiculairement a la 
direction de Taxe principal 19 de Timage suivant lequel se produit Talteration de 
couleur. La reunion de Tensemble des pixels de chaque bande 18 represente 

10 Tensemble des pixels de Timage 12. La division de Timage 12 en bandes de pixels 
18 est operee suivant des bandes 18 paralleles adjacentes; selon un premier mode 
de realisation, ces bandes paralleles adjacentes n'ont aucun pixel commun entre 
elles, ou bien selon un second mode de realisation, la division de Timage est op6ree 
suivant des bandes paralleles se recouvrant partiellement deux a deux suivant au 

15 moins une rangee de pixels; un troisieme mode de realisation permet de diviser 
Timage a la fois en bandes adjacentes ou recouvrantes. Le mode de realisation 
choisi parmi ces trois variantes depend de Valgorithme utilis6 pour d6placer, dans 
Timage 12, la bande mobile 1 8 de largeur 21 . Dans un mode de realisation 
preferentiel, les bandes 18 sont de meme largeur 21. Mais Talgorithme mettant en 

20 oeuvre le procede de Tinvention permet egalement de choisir des bandes 1 8 de 
largeurs 21 inegales. 

A partir de la division de Timage num6rique 12 en bandes de pixels 
1 8, le procede de Tinvention permet de calculer automatiquement, pour chaque 
bande de pixels 18, et selon la figure 6, pour chacun des plans des canaux colores 

25 40, 50, 60 foimant Timage 12 (par exemple les canaux rouge, vert et bleu), une 

distribution des densites optiques de chaque pixel de ladite bande 18. Un exemple 
de distribution des densites optiques est represents a la figure 4. 

La d6tennination des variations spatiales de Talteration des couleurs 
dans Tensemble de Timage 12, et pour chacun des canaux de couleur, par exemple 

30 rouge, vert, bleu, consiste a separer ou filtrer la contribution du contenu de Timage, 
de la contribution de Taltdration des couleurs. Ceci en tenant compte de la variation 
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de la distribution des valeurs codecs de densite optique, comme represente dans la 
figure 4. 

La figure 5 represente, pour un canal de couleur de l'image, une 
representation spatiale de la distribution reelle des densites optiques dans 
5 l'extension Xe de ladite image. L'extension Xe represente la dimension de l'image 
suivant laquelle se produit Alteration ou degradation de couleur (axe principal 19). 
Les lignes 41 et 45 correspondent aux bords de Timage. La distribution reelle des 
densites optiques est representee par exemple par les courbes 30, 3 1, 32, 33, 34. 
Ces courbes 30, 3 1, 32, 33, 34 correspondent a un mixage du profil de degradation 

10 des couleurs et du contenu de Timage. Le proc<§de de 1'invention permet de separer 
automatiquement les signaux lies au contenu de l'image d'une part, et au profil de 
degradation des couleurs d'autre part. 

Selon la figure 5, le proced6 de 1'invention permet de calculer 
automatiquement, pour chacune des distributions de densites optiques des pixels 

15 de la bande de pixels 44, une valeur de reference brute superieure de densite 

optique 36, et une valeur de reference brute inferieure de densite optique 37. Selon 
la figure 5, la distribution des densites optiques est representee par une famille de 
courbes 30, 31, 32, 33, 34. Les courbes 30 et 34 represented par exemple les 
courbes de densite optique respectivement du maximum et du minimum de la 

20 distribution. Dans un mode de realisation prdfere, ralgorithme du procede de 
rinvention ^limine automatiquement ces deux extremes de la distribution (les 
courbes 30 et 34) pour operer les calculs. Les valeurs de reference brute superieure 
36 et inferieure 37 sont calculees a partir des courbes 3 1 et 32 qui correspondent 
par exemple respectivement aux courbes de densite optique de la famille de 

25 courbes plac£es dans la distribution, au voisinage respectivement des courbes 
extremes 30 et 34. 

Toutefois, dans un autre mode de realisation, les valeurs de 
reference brutes superieures de densite optique 39, et inferieures 38, sont calculees 
h partir des courbes de densite optique du maximum 30 et du minimum 34. 
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A partir des valeurs de reference brutes superieure 36 et inferieure 
37 de densite optique, le procede de l'invention permet de calculer des valeurs de 
reference filtr6es superieure 46 et inferieure 47 de densite optique. 

Le procede de l'invention, a partir des valeurs de reference filtrees 
46, 47, determine automatiquement des courbes de reference de densite optique: la 
courbe de reference sup6rieure de densite optique 42, et la courbe de reference 
inferieure de densite optique 43. Les courbes de reference de densite optique 42, 43 
sont representatives du profil de degradation de couleur de 1'image. Ces courbes de 
reference 42, 43 sont obtenues par filtrage, en tenant compte de la quatrieme 
hypothese liee a la mise en ceuvre de l'invention. La quatrieme hypothese est basee 
sur un profil de degradation de couleur de Hmage qui varie de maniere continue, 
Hsse, et monotone. Dans un mode de realisation prefere, reoperation de filtrage, 
dont le principe est represents sur la figure 9, est operee de maniere a appliquer 
cette quatrieme hypothese, eliminer les points extremes de la distribution 
correspondent, apres filtrage, aux courbes extremes de la distribution Co et Cioo, et 
separer le contenu de 1'image du profil de degradation de couleur de ladite image. 
Les courbes C 0 et Cioo sont les courbes de r6ference inferieure et superieure de 
densite optique qui represented les courbes extremes de la distribution. 

Selon la figure 9 et dans le mode de realisation prefere de 
l'invention, l'operation de filtrage est operee en utilisant la methode dite du suivi 
des maximum. Cette methode de filtrage revient, a partir d'une courbe 3 1 
representative d'une distribution, a obtenir la courbe resultante qui represente le 
minima des deux courbes 35C et 36C. Ceci permet d'eliminer les creux de la 
courbe initiale 31 et ainsi de respecter l'hypothese de monotonia Selon la figure 6, 
pour chacun des canaux de couleur 40, 50, 60 formant 1'image, les valeurs de 
reference filtrees superieures de densite optique 46, 56, 66, et les valeurs de 
reference filtrees inferieures de densite optique 47, 57, 67, sont calculees a partir 
de rutilisation d'un filtrage de 1'image 12 par la methode du suivi des maximum 
dans chaque bande de pixels 18. 

Selon les figures 6 et 10, le procede de l'invention permet ainsi 
d'obtenir, apres filtrage, des courbes corrigees lissees 42, 43, 52, 53, 62, 63. 
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Suivant un mode prefere de realisation, les courbes C 0 et C 10 o, correspondant aux 
extremes de la distribution (0% et 100%), ont ete exclues. La figure 6 represente 
les courbes de reference superieures de densit6 optique 42, 52, 62 et 
respectivement inferieures 43, 53, 63. Les dites courbes de reference etant 
representees dans chaque plan 40, 50, 60 correspondant aux canaux de couleur. Par 
analogie avec la representation graphique de la figure 3, l'extension de Timage 
suivant laquelle se produit la variation de Talteration de couleur est representee en 
abscisse sur la figure 6 par la distance comprise entre les lignes 41 et 45 pour le 
canal de couleur 40, entre les lignes 45 et 55 pour le canal de couleur 50, et entre 
les lignes 55 et 65 pour le canal de couleur 60. L'axe des ordonnees de la figure 6 
represente les valeurs codees des densites optiques V c (DO). La couche d'emulsion 
photographique de Timage 3 correspondant a la couleur rouge est representee en 
abscisse par exemple par le plan 40 de Timage; la couche d'emulsion 
photographique correspondant a la couleur verte est representee en abscisse par 
exemple par le plan 50 de Timage; La couche d'emulsion photographique 
correspondant a la couleur bleue est representee en abscisse par exemple par le 
plan 60 de Timage. Selon la figure 6, le procede de Tinvention permet de 
determiner automatiquement, pour chacun des canaux de couleur 40, 50, 60 
formant Timage, et a partir de la valenr maximale de densite optique 46M, 56M, 
66M de la courbe de reference superieure de densite optique, la bande de pixels de 
moindre degradation 44M, 54M, 64M. Les valeurs de reference filtrees superieures 
de densite optique des bandes de pixels de moindre degradation 44M, 54M, 64M 
sont respectivement 46M, 56M, 66M. Les valeurs de reference filtrees inferieures 
de densite optique des bandes de pixels de moindre degradation 44M, 54M, 64M 
sont respectivement 47M, 57M, 67M. Cette bande de pixels de moindre 
degradation correspond a la bande de pixels de Timage ayant subi une alteration ou 
degradation minimale de couleur. Selon un mode parti oilier de realisation du 
procede de Tinvention, la bande de moindre degradation est la bande de pixels 
comprenant une seule rangee de pixels, et passant par le centre 20 de Timage 12. 
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La figure 7 represente la fonction de restauration des densites 
optiques. La transformation lineaire est operee suivant une equation de 
transformation Eq(2) du type: Dr = a Dp + b; 

Equation dans iaqueile a et b sont independants pour chaque canal colore rouge (a r , 
b r ), vert (a v , b v ), bleu (a^ bb), et dans Iaqueile a et b sont determines par 
ralgorithme permettant de mettre en oeuvre le procede de Tinvention. Cette 
seconde equation Eq(2) represente la transformation inverse de celle de la premiere 
equation Eq(l). L'application simultanee de Tequation Eq(2) aux trois canaux 
colores s'exprime sous la forme: 

DRr = a r DFr + b r 

Dry 1=5 a v Dpv + b v 
Dju> = a b DFb + b b 
soit, sous la forme matricielle: 



Drt 








Dju, 
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D Fr 
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0 




D Fy 
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b v 


0 


0 


a b 




Da 




b b 



Dans tin mode de realisation particulier de Tinvention, les termes 
non diagonaux de la matrice a trois lignes et trois colonnes, peuvent etre non mils 
pour tenir compte de l'effet du recouvrement spectral des absorptions, conune 
decrit par exemple dans le brevet US 5 796 874 prScite. 

Selon les figures 6, 7, et 10, le procSde de l'invention pennet ensuite 
d'appliquer a chaque bande de pixels 44, 54, 64 de Timage numerique autre que la 
bande de pixels de moindre degradation 44M, 54M, 64M, et pour chacun des 
canaux de couleur formant I'image, une transformation lineaire permettant de 
transformer les valeurs filtrees de densite optique 46, 56, 66 placees sur la courbe 
de reference superieure de densite optique 42, 52, 62 et les valeurs filtrees de 
densite optique 47, 57, 67 placees sur la courbe de reference inferieure de densite 
optique 43, 53, 63 de ladite chaque bande de pixels 44, 54, 64, respectivement en 
valeurs sup6rieures 46M, 56M, 66M et inferieures 47M, 57M, 67M de densite 
optique de la bande de pixels de moindre degradation 44M, 54M, 64M. 
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Le precede de Tinvention permet done d'obtenir automatiquemeat 
une image numerique 12 dont les couleurs sont restaurees, e'est-a-dire dont les 
couleurs ont ete transformees pour corresponds a celles de Timage argentique 
d'origine. Uoperateur peut par exemple, apres avoir applique le procede de 
5 Tinvention sur une image alteree, afficher sur l'ecran 1 1 l'image numerique 1 2 dont 
les couleurs sont restaurees; il peut egalement restituer ensuite cette image 
restauree, par exemple sur du papier en utilisant une imprimante laser ou a jet 
d'encre, ou bien sur un film photographique. 

L'invention a ete decrite en faisant reference en particulier a des 
10 modes de realisation preferes de celle-ci; mais il est evident que des variantes et 
modifications peuvent etre apportees a Tinterieur de la port6e des revendications. 
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REVENDICAHONS 
Proced6 automatique de restauration de la couleur d'une image argentique 
(3) dont la couleur est alteree par rapport a la couleur d'origine de ladite 
image argentique (3), l'alteration de la couleur se produisant non 
umformement, spatialement dans le plan de l'image, et suivant une 
direction parallele a un axe principal (19) de l'image, ledit procede 
comprenant les etapes suivantes: 

a) numeriser 1'image argentique alteree en uulisant une echelle de 
mesure des densites optiques; 

b) filtrer l'ensemble des pixels de l'image numerisee; 

c) diviser l'image numerique (12) en bandes de pixels (18) 
disposees perpendiculairement a la direction de l'axe principal (19) de 
l'image suivant lequel se produit l'alteration de couleur, la reunion de 
l'ensemble des pixels de chaque bande representant l'ensemble des pixels 
de l'image numerique; 

d) calculer, pour chaque bande de pixels, et pour chacun des canaux 
de couleur (40), (50), (60) formant l'image, une distribution des densites 
optiques des pixels formant ladite bande; 

e) calculer, pour chacune desdites distributions de l'etape d), une 
valeur de reference brute superieure de densite optique (36), et une valeur 
de ref6rence brute inferieure de densite optique (37); 

f) calculer, a partir de l'ensemble des valeurs de reference brutes 
superieures et inferieures de l'etape e), respectivement des valeurs de 
reference filtrees superieures de densite optique (46) et des valeurs de 
reference filtrees inferieures de densite optique (47) correspondantes; 

g) determiner, a partir des valeurs de reference filtrees sup&ieures 
et inferieures de densite optique de l'etape f), respectivement des courbes d< 
reference superieures de densite optique (42) et des courbes de reference 
inferieures de densite optique (43), lesdites courbes dereference etant 
representatives du profil de degradation de couleur de l'image; 
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h) determiner, pour chacun des canaux de couleur fbrraant I'image, 
et a partir de la valeur maximale de densite optique (46M), (56M), (66M) 
de la courbe de reference superieure de densite optique (42), la bande de 
pixels demoindre d6gradation (44M), (54M), (64M); 

i) appliquer a chaque bande de pixels (44), (54), (64) de I'image 
numerique autre que la bande de pixels de moindre d6gradation (44M), 
(54M), (64M), et pour chacun des canaux de couleur formant I'image, une 
transformation lineaire permettant de transformer les valeurs filtrees de 
densite optique (46), (56), (66), (47), (57), (67) placees sur la courbe de 
reference superieure (42), (52), (62) et sur la courbe de reference inferieure 
(43), (53), (63) de densite optique de ladite chaque bande de pixels (44), 
(54), (64), respectivement en valeurs superieure (46M), (56M), (66M) et 
inferieure (47M), (57M), (67M) de densite optique de la bande de pixels de 
moindre degradation. 

Proced6 selon la revendication 1, dans lequel le filtrage initial de I'image 
numerique alteree est opere en unlisant des operateurs de morphologie 
mathematique du type ouverture d'ordre (n) suivie de fermeture d'ordre (n), 
(n) &ant un nombre entier. 

Precede selon 1'une quelconque des revendications 1 ou 2, dans lequel les 
valeurs de reference filtrees superieures de densite optique (46), (56), (66), 
et les valeurs de reference filtrees inferieures de densite optique (47), (57), 
(67), sont calculees a partir de l'utilisation d'un filtrage de I'image (12) par 
la methode du suivi des maximum dans chaque bande de pixels (1 8). 

Proc6de selon l'une quelconque des revendications 1 a 3, dans lequel la 
division de I'image numerique (12) en bandes de pixels (1 8) est operee 
suivant des bandes paralleles adjacentes n'ayant aucun pixel commun entre 
elles. 
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Procede selon l'une quelconque des revendications 1 a 3, dans lequel la 
division de 1'image numerique (12) en bandes de pixels (18) est operee 
suivant des bandes paralleles se recouvrant partiellement deux a deux 
suivant au moins une rangee de pixels. 

Proced6 selon l'une quelconque des revendications precedentes, dans lequel 
les valeurs de reference brutes superieures (36) et inferieures (37) de 
densite optique correspondent respectivement aux courbes de la 
distribution de densite optique (31), (32) placees au voisinage des courbes 
de densite optique maximum (30) et minimum (34) de ladite distribution. 

Procede selon l'une quelconque des revendications 1 a 5, dans lequel les 
valeurs de reference brutes superieures (39) et inferieures (38) de densite 
optique correspondent respectivement aux courbes de la distribution de 
densite optique maximum (30) et niinimum (34) de ladite distribution. 

Procede selon l'une quelconque des revendications precedentes, dans lequel 
la bande de moindre degradation (44M), (54M), (64M) est une bande de 
pixels de l'image comprenant le pixel de valeur maximale de densite 
optique (46M), (S6M), (66M) de la courbe de reference superieure de 
densite optique (42), (52),(62). 

Procede selon l'une quelconque des revendications precedentes, dans lequel 
la bande de moindre degradation est la bande de pixels comprenant une 
seule rangee de pixels, ladite bande passant par le centre (20) de l'image 
(12). 

Procede selon l'une quelconque des revendications precedentes, dans lequel 
les plans des canaux de couleur formant l'image sont le rouge, le vert et le 
bleu. 
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Procede selon Tune quelconque des revendications 1 a 3, dans lequel la 
division de l'image numerique (] 2) en bandes de pixels (] 8) est operee 
suivant des bandes paralleles se recouvrant partiellement deux a deux 
suivant au moins une rangee de pixels. 

Procede selon I'une quelconque des revendications precedentes, dans lequel 
les valeurs de reference brutes superieures (36) et inferieures (37) de densite 
optique correspondent respectivement aux courbes de la distribution de 
densite optique (3 1), (32) placees au voisinage des courbes de densite 
optique maximum (30) et minimum (34) de Iadite distribution. 

Procede selon Tune quelconque des revendications 1 a 5, dans lequel les 
valeurs de reference brutes superieures (39) et inferieures (38) de densite 
optique correspondent respectivement aux courbes de la distribution de 
densite optique maximum (30) et minimum (34) de ladite distribution. 

Procede selon I'une quelconque des revendications precedentes, dans lequel 
la bande de moindre degradation (44M), (54M), (64M) est une bande de 
pixels de l'image comprenant le pixel de valeur maximale de densite optique 
(46M), (56M), (66M) de la courbe de reference superieure de densite 
optique (42), (52),(62). 

Procede selon I'une quelconque des revendications 1 a 7, dans lequel la 
bande de moindre degradation est la bande de pixels comprenant une seule 
rangee de pixels, ladite bande passant par le centre (20) de J'image (12). 

Procede selon I'une quelconque des revendications precedentes, dans lequel 
les plans des canaux de couleur formant l'image sont le rouge, le vert et le 
bleu. 
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